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Abstract. Partition functions for dimers on closed oriented surfaces are known 
to be alternating sums of Pfaffians of Kasteleyn matrices. In this paper, we 
obtain the formula for the coefficients in terms of discrete spin structures. 
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Introduction 

Dimcr models on graphs have a long history in statistical mechanics [HI I15j . 
States in dimer models are perfect matchings between vertices of the graph where 
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only adjacent vertices are matched. The probability of a state is determined by 
assigning weights to edges. 

Dimer models also have many interesting mathematical aspects involving com- 
binatorics, probability theory [111 real algebraic geometry jlOl etc... One 
of the remarkable facts about dimer models is that the partition function can be 
written as a linear combination of 2 2g Pfaffians of N x N matrices, where N is the 
number of vertices in the graph and g the genus of a surface where the graph can 
be embedded. 

The matrices in the Pfaffian formula for the dimer partition function are called 
Kasteleyn matrices. They involve certain orientations of edges of the graph known 
as Kasteleyn orientations. Two Kasteleyn orientations are called equivalent if one 
can be obtained from the other by a sequence of moves reversing orientations of all 
edges adjacent to a vertex. 

The number of non-equivalent Kasteleyn orientations of a surface graph of genus 
g is 2 2g and is equal to the number of non-equivalent spin structures on the surface. 
The Pfaffian formula expresses the partition function of the dimer model as an 
alternating sum of Pfaffians of Kasteleyn operators, one for each equivalence class 
of Kasteleyn orientations. This formula was proved in for the torus, and it was 
stated in that for other surfaces, the partition function of the dimer model is 
equal to the sum of 2 2g Pfaffians. The formal combinatorial proof of this fact and 
the exact description of coefficients for all oriented surfaces first appeared in ^5] 
and ^Hj ( see a l so E9)- A combinatorial proof of such formula for non-orientable 
surfaces can also be found in 

The partition function of free fermions on a Ricmann surface of genus g is also a 
linear combination of 2 2g Pfaffians of Dirac operators. Each term in this sum cor- 
responds to a spin structure pQ| . Assuming that dimer models are discretizations of 
free fermions on Riemann surfaces, one should expect a relation between Kasteleyn 
orientations and spin structures and between the Kasteleyn operator for a given 
Kasteleyn orientation and the Dirac operator in the corresponding spinor bundle. 
Numerical evidence relating the critical dimer model on a square and triangular 
lattices in the thermodynamical limit with Dirac operators can be found in |2] for 
.9 = 2. 

An explicit construction relating a spin structure on a surface with a Kasteleyn 
orientation on a graph with dimer configuration was suggested in Furthermore, 
for bipartite graphs with critical weights, the Kasteleyn operator can be naturally 
identified with a discrete version of the Dirac operator [5] . This gives an interesting 
relation between dimer models and the theory of discrete meromorphic functions 
d 

In this paper we investigate further the relation between Kasteleyn orientations 
and spin structures, and use this relation to give a geometric proof of the Pfaffian 
formula for closed surfaces. Below is a brief summary of our main results. 

Recall some basic notions. A dimer configuration on a graph T is a perfect 
matching on vertices where matched vertices are connected by edges. Given two 
such configurations D and D', set A(D, D') = (DUD 1 ) \(Dnfl'). A surface graph 
is a graph T embedded into a surface £ as the 1-squeletton of a CW-decomposition 
of S. A Kasteleyn orientation of a surface graph is an orientation of edges of the 
graph, such that the product or relative orientations of boundary edges of each face 
is negative (see Section EJ. 
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One of our results is that any dimer configuration Dona surface graph r C £ 
induces an isomorphism of affine -ff^E; Z2)-spaces 

(1) i, D : (DC(T)/~) —> 0(^(2; Z 2 ), •) , [#] — » <?£ 

from the set of equivalence classes of Kasteleyn orientations on Y C E onto the set 
of quadratic forms on (Hi (E; Z2), •), where ■ denotes the intersection form on E. 
Furthermore, ipo = V'D' if and only if Z? and D' are equivalent dimer configurations 
(that is, A(D, D') is zero in fli(E;Z 2 )). 

Since the affine space of spin structures on S is canonically isomorphic to the 
affine space of such quadratic forms, this establishes an isomorphism between equiv- 
alence classes of Kasteleyn orientations and spin structures. 

This correspondence implies easily the following identity. Let Do be a fixed 
dimer configuration on a graph Y. Realize T as a surface graph T C E of genus 
g. Let K be any Kasteleyn orientation on Y C E, and let A K be the associated 
Kasteleyn matrix. Then, 

(2) Yi(A K ) = e K (D a ) (-l) qE ° ia) Za(D ), 

QeHi(S;Z 2 ) 

where q^ g is the quadratic form associated to K and Do via 0J, e (Do) is some 
sign depending on K and Do, and 

Z a (D Q ) = J2 W ( D )' 

D 

the sum being on all dimer configurations D such that A(Dq,D) is equal to a in 
fli(£;Za). 

It follows that the partition function of a dimer model on F is given by 

(3) Z= i£Ar%5 o ) £ *(Do)Pf(^), 

[K] 

were the sum is taken over the 2 2g equivalence classes of Kasteleyn orientations on 
F C E, and Arf(q^ ) = ±1 denotes the Arf invariant of the quadratic form q^ Q . 
Note that the sign Ari(q^ g )e K (D ) does not depend on D - 

The paper is organized as follows. In Section ^ we introduce the dimer model 
on a graph Y and define composition cycles. Section |3 deals with dimers on surface 
graphs r C E and the definition of an equivalence relation for dimer configurations 
on surface graphs. In Section |3 we recall the definition of a Kasteleyn orientation. 
We then show that a surface graph F C E admits such an orientation if and only 
if the number of vertices of Y is even. We prove that, in such a case, the set 
of equivalence classes of Kasteleyn orientations on Y C E is an affine i/ 1 (E;Z 2 )- 
space. Finally, we give an algorithmic procedure for the construction of the 2 2g 
non-equivalent Kasteleyn orientations on a given surface graph T C E of genus 
g. The core of the paper lies in Section where we establish the correspondence 
stated above. This result is used in Section |S] to obtain equations © and 
We also give a formula for the local correlation functions of a dimer model. In 
the appendix, we collect formulae expressing dimer models in terms of Grassman 
integrals. 
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1.1. Dimer configurations and composition cycles on graphs. Let T be a 

finite connected graph. A perfect matching on T is a choice of edges of T such that 
each vertex of T is adjacent to exactly one of these edges. In statistical mechanics, 
a perfect matching on T is also known as a dimer configuration on T. The edges of 
the perfect matching are called dimer s. An example of a dimer configuration on a 
graph is given in Figure ^ 



In order to have a perfect matching, a graph T clearly needs to have an even 
number of vertices. However, there are connected graphs with an even number 
of vertices but no perfect matching. We refer to ^3] for combinatorial aspects of 
matchings. Throughout the paper and unless otherwise stated, we will only consider 
finite graphs which admit perfect matchings. In particular, all the graphs will have 
an even number of vertices. 

Given two dimer configurations D and D' on a graph T, consider the subgraph of 
r given by the symmetric difference (DU D')\{D DD') . The connected components 
of this subgraph are called (D, D')- composition cycles or simply composition cycles. 
Clearly, each composition cycle is a simple closed curve of even length. This is 
illustrated in Figure |3 where the two dimer configurations are shown in black and 
traced lines. 

1.2. Edge weight system. Let D(T) denote the set of dimer configurations on a 
graph r. A weight system on D(T) is a positive real- valued function on this set. A 
weight system w defines a probability distribution on all dimer configurations: 



1. The dimer model 




Figure 1 . A dimer configuration on a graph. 



Prob(D) 



w(D) 
Z(T;w) 



where 



Z(T-w) 



D 



DIMERS ON SURFACE GRAPHS AND SPIN STRUCTURES. I 



5 




Figure 2. An example of composition cycles: dimers from D are 
in solid, and dimers from D' are in traced lines. On this example, 
there is one (D, Decomposition cycle of length 2, one of length 4, 
and one of length 6. 



is the partition function. This probabilistic measure is the Gibbs measure for the 
dimer model on the graph T with the weight system w. 

We shall focus on a particular type of weight system called edge weight system. 
Assign to each edge e of T a positive real number w(e), called the weight of the 
edge e. The associated edge weight system on D(T) is given by 

w(D) = J] «>(e), 

eeD 

where the product is over all edges occupied by dimers of D. In statistical mechan- 
ics, these weights are called Boltzmann weights. Their physical meaning is 

w(e) = exp(^- -y-J, 

where E(e) is the energy of dimer occupying the edge e and T is the temperature. 

1.3. Local correlation functions. Let e be an edge of T. The characteristic 
function of e is the function a e on D(T) given by 



O-e(D) 



1 if eeD; 
otherwise. 



The expectation values of products of characteristic functions are called local 
correlation functions, or dimer-dimer correlation functions: 

Z(ei,...,e k ;T;w) 

— — ' 

where 

k 

Z(e 1 ,...,e k ;T-w)=Y / I[^)]l^(D)= ]T w(D). 

D e£D i=l D3e!,...,e k 

Note that < a ei ■ ■ -a ek >= if some edges ei ^ ej share a common vertex. Note 
also that a e ■ a e — a e . Therefore, it may be assumed that ej ^ ej for i ^ j. 
For any dimer configuration D, 

w(D) 



/ TT rr \ W< " L 

eeD 
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so we can reconstruct the weight system if we know all local correlation functions. 
In this sense, local correlation functions carry all the information about the Gibbs 
measure. 

2. DlMERS ON SURFACE GRAPHS 

2.1. Surface graphs. Let E be a connected oriented closed surface. By a surface 
graph, we mean a graph T embedded in E as the 1-squeleton of a cellular decompo- 
sition X of E. We shall assume throughout the paper that the surface E is endowed 
with the counter-clockwise orientation. 

Any hnite connected graph can be realized as a surface graph. Indeed, such a 
graph r always embeds in a closed oriented surface of genus <?, for g sufficiently large. 
If the genus is minimal, one easily checks that T induces a cellular decomposition 
of the surface E. 

We will be interested mostly in large graphs. In this paper we will focus on 
graphs embedded into a surface of fixed genus. 

2.2. Equivalent dimer configurations. A dimer configuration on a surface graph 
r C E is simply a dimer configuration on the graph T. Such a dimer configuration 
can be regarded as a 1-chain in the cellular chain complex of A" with Z2-coefficients: 

c D = ^eeC7 1 (A;Z 2 ). 

egfl 

By definition, dcu = St,er v e Co(A;Z 2 ), the sum being on all vertices v of T. 
Therefore, given any pair of dimer configurations D and D' on T C E, Cd + Cd< is 
a 1-cycle: 

d(c D +c D ,)= 8c D + dc D , = + ") = »£ 0)(A; Z 2 ). 

This 1-cycle is nothing but the union of all (D, D')-composition cycles. Let A(D, D') 
denote its homology class in i?i(X;Z 2 ) = iii(E;Z 2 ). We shall say that two dimer 
configurations D and D' are equivalent if A(D, D') = in i?i(E; Z 2 ). 

Note that these concepts make perfect sense when T is the 1-squeleton of any 
CW-complex, not necessarily the cellular decomposition of an oriented closed sur- 
face. 

3. Kasteleyn orientations on surface graphs 

Let r C E be a surface graph. The counter-clockwise orientation of E induces 
an orientation on each 2-cell, or face of X. An orientation K of the edges of V is 
called a Kasteleyn orientation if for each face / of X, 

(4) n £ f ( e ) = -L 

where the product is taken over all boundary edges of /, and 

K( . j 1 if e is oriented by K as the oriented boundary of the face /; 
£ f ye) = s 

1—1 otherwise. 

This is illustrated in Figure 03 

Define the operation of orientation changing at a vertex as the one which flips 
the orientation of all the edges adjacent to this vertex, as illustrated in Figure 0] 



DIMERS ON SURFACE GRAPHS AND SPIN STRUCTURES. I 



7 




Figure 3. A Kasteleyn orientation on the boundary edges of a face. 

It is clear that such an operation brings a Kasteleyn orientation to a Kasteleyn 
orientation. Let us say that two Kasteleyn orientations are equivalent if they are 
obtained one from the other by a sequence of orientation changes at vertices. 




Figure 4. Orientation change at a vertex. 



3.1. Existence of a Kasteleyn orientation. 

Theorem 3.1. There exists a Kasteleyn orientation on a surface graph T C £ if 
and only if the number of vertices of T is even. 

Proof. Let uj be any orientation of the edges of T, and let c w £ C 2 (X; Z2) be defined 
by the equation 

(_!)«*(/) = _ JJ £ - (e ). 
eedf 

Note that u> is Kasteleyn if and only if c w = 0. Let V, E and F denote the 
number of vertices, edges and faces in X, respectively. Since E is closed, its Euler 
characteristic is even, and we get the equality mod 2 

- X (E) = V + E + F = V+J2 C M)- 

Here, each edge e contributes to the number c w (/) where / is the face whose oriented 
boundary contains e with the orientation opposite to uj. Therefore, the number of 
faces / such that c UJ {f) — 1 has the parity of V. Hence, if V is even, then c w (/) = 1 
for an even number of faces, so c w is a 2-coboundary. In other words, c w = 5a for 
some a £ C 1 (X; Z2). Let K be the orientation of the edges of L which agrees with 
lu on e if and only if <r(e) = 0. Clearly, ck = 0, so K is a Kasteleyn orientation. 
Conversely, let us assume that there is a Kasteleyn orientation K . This means that 
ck (/) = 1 for none of the faces. By the argument above, V is even. □ 

A more constructive proof of this result will be given in Section 13.31 
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3.2. Uniqueness of Kasteleyn orientations. Let V be a vector space. Recall 
that an affine V -space is a set S endowed with a map S x S — > V, (a,b) ^ a — b 
such that: 

(i) for every a, b and c in S, we have (a — b) + (b — c) = a — c ; 

(ii) for every b in S, the map S — > V given by a i— > a — 6 is a bijection. 

In other words, an affine V^-space is a V-torsor: it is a set endowed with a freely 
transitive action of the abelian group V. 

Theorem 3.2. Let T C £ &e a surface graph with an even number of vertices. 
Then, the set of equivalence classes of Kasteleyn orientations onTcS is an affine 
-ff Z 2 )-sj?ace. 

Corollary 3.3. There are exactly 2 2g equivalence classes of Kasteleyn orientations 
on F C E, where g denotes the genus of E. □ 

Proof of the theorem. Let 3C(T) denote the set of Kasteleyn orientations on T C £. 
By Theorem mi it is non-empty. Consider the map d: %{T) x X(T) -> C X (X; Z 2 ) 
given by $k.k' (e) = if if and if' agree on the edge e, and $k,K' (e) = 1 otherwise. 
Since if and if' are Kasteleyn orientations, 

( _l)*^(»/3 = JJ (_i)*^W = J| ^(e) . ft £ f (e) = ( -i)(_i) = i 

for any face /. Therefore, 5*&K,K'(f) = $K,K'(df) = 0, that is, $k,K' is a 1-cocycle. 

Thus, we get a map DC(r) x 0C(r) ^ H 1 ^^). Note that ^.Jf +i?if',K" =$k.k>> 
for any if, if', if" e 3C(r). Also, one easily checks that = in iJ^EjZa) 

if and only if there is a sequence of vertices such that if is obtained from if' by 
reversing the orientation around all these vertices, i.e., if and only if if ~ if'. It 
follows that we have a map 

(oc(r)/ ~) x (ac(r)/ ~) — h 1 ^-, z 2 ), ([if], \k'\) » [if] - \k'\ ■.= [t KiK ,] 

such that for any K' in 3C(r), the map (3C(r)/ ~) -> ii 1 (S;Z 2 ) given by [if] ^ 
[if] — [if] is injective. Finally, let us check that this map is onto. Fix a class in 
ii 1 (£;Z2), and represent it by some 1-cycle a 6 Z 1 (X;Z 2 ). Let if be the same 
orientation as if' whenever <r(e) = 0, and the opposite when <r(e) = 1. Obviously, 
"&k,k' = <y. Furthermore, if is Kasteleyn since if' is and 5a = 0. Indeed, given a 
face /, 

0=(5a)(f)=a(df)=ti K , K ,(df), 
so 1 = Ueedf(~ 1 f K - K ' {e) = (- 1 ) Ueedf £ f( e )- This concludes the proof. □ 

3.3. How to construct Kasteleyn orientations. Given a surface graph r C S 
with an even number of vertices, we know that there are exactly 2 2s non-equivalent 
Kasteleyn orientations on T. However, the proof given above is not really con- 
structive. For this reason, we now give an algorithm for the construction of these 
Kasteleyn orientations. The successive steps of the algorithm are illustrated in 
Figure 

(0) Let r C S be a surface graph with an even number of vertices, and let g 
denote the genus of E. 

(1) Consider a system a = a\ U- • • Ua 2g of simple closed curves on T such that 
E cut along a is a 2-disc E'. (Note that such curves exist since T induces 
a cellular decomposition of E.) 
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(2) The surface graph Y C E induces a graph V C E' which is the 1-squeleton 
of a cellular decomposition of the 2-sphere S 2 . Fix a spanning tree T of the 
graph dual to the surface graph f C S 2 , rooted at the vertex corresponding 
to the face S 2 \E'. 

(3) Orient the edges of V which do not intersect T respecting the following 
condition: whenever two edges of V in <9E' are identified in E, their ori- 
entations agree. Now, every edge in 9E' is oriented, except one. Let us 
denote it by e*. 

(4) Orient edges of T' such that the Kasteleyn condition holds for the faces 
corresponding to the leaves of T. Moving down the tree (from the leaves 
to the root), orient each crossing edge such that the Kasteleyn condition 
holds for all faces left behind. 

This gives a Kasteleyn orientation of f c S'. By the condition in step 3, it induces 
a Kasteleyn orientation on Y C E, provided that the orientation of the last edge 
e» satisfies this condition. It turns out to be the case if and only if the number 
of vertices of T is even. (This is an easy consequence of Theorem 13.11 ) Therefore, 
we have constructed a Kasteleyn orientation K on Y C E. To obtain the 2 2g 
non-equivalent ones, proceed as follows. 

(5) Consider a family of simple closed curves j3\,...,j32 g on E avoiding the 
vertices of Y, and forming a basis of Zfi(E; Z2). 

(6) Consider the Kasteleyn orientation K, and some subset / C {1, . . . ,2g}. 
For all i € I, change the orientation of all the edges in Y that intersect 

The resulting orientation Kj is clearly Kasteleyn, as df ■ (3i is even for every face / 
and index i. Furthermore, one easily checks that Ki and Kj are non-equivalent if 
I J. Hence, we have constructed the 2 2g non-equivalent Kasteleyn orientations 
on T C E. 




Figure 5. An example of the explicit construction of a Kasteleyn orientation. 

4. Kasteleyn orientations as discrete spin structures 

We saw in Corollary 13.31 that there are exactly 2 2ff non-equivalent Kasteleyn 
orientations on a surface graph T C E, where g denotes the genus of E. It is known 
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that this is also the number of non-equivalent spin structures on E. This relation 
between the number of Kasteleyn orientations and the number of spin structures 
is not accidental. Kasteleyn orientations of surface graphs can be regarded as 
discrete versions of spin structures. This statement will be made precise in the 
present section. (See in particular Corollary 14. 41 ) 

4.1. The quadratic form associated to a Kasteleyn orientation. Let V be 

a finite dimensional vector space over the field Z2, and let ip: V x V — > Z2 be a 
fixed bilinear form. Recall that a function q : V — > Z2 is a quadratic form on (V, ip) 
if 

q(x + y) = q{x) + q(y) + <p(x, y) 
for all x,y G V. Note that the difference (that is, the sum) of two quadratic forms 
on (V, ip) is a linear form on V. Therefore, one easily checks that the set Q(V, ip) of 
quadratic forms on (V, ip) is an affine V*-spa.ce, where V* denotes the dual of V . 

Fix a Kasteleyn orientation if on a surface graph r G E. Given an oriented 
simple closed curve C on T, set 

eeC 

where £ci e ) ls equal to +1 (resp. —1) if the orientations on the edge e given by 
C and K agree (resp. do not agree). For a fixed dimer configuration D on L, let 
£d(C) denote the number of vertices v in C whose adjacent dimer of D sticks out 
to the left of C in E. 

Theorem 4.1. Given a class a £ 7?i(E;Z2) ; represent it by oriented simple closed 
curves C\,..., C m in V . If K is a Kasteleyn orientation ofiTcS, then the function 
q% : iii(E; Z 2 ) — > Z 2 given by 

m 

(_i)<£(«) = (-!)£.<, c -^;Q(- e ^(a))(-i) £D(Ci) 

1=1 

is a well-defined quadratic form on (iii(E; Z2), ■), where ■ denotes the intersection 
form. 

We postpone the proof of this result to the next subsections. Let us first inves- 
tigate some of its consequences. 

Proposition 4.2. (?) Let D be a fixed dimer configuration on T. If K and K' are 

two Kasteleyn orientations on T C E, then q^ — q^j maps to [K] — [K'] via the 
canonical isomorphism Hom(7?i(E; Z2); Z2) = i/ 1 (E;Z2). 

(ii) Let K be a fixed Kasteleyn orientations on T C E. If D and D' are two 
dimer configurations on T, then q^ — q^, G Hom(i?i(E; Z2); Z2) is given by a t—> 
a ■ A(D,D'). 

Proof. Let C be a simple closed curve in T representing a class a in i?i(E; Z2). By 
definition, 

(_!)(<£-<£')(«) = Y[ 4(e) 11 e%\e) = J](-l)^.-'( e ) = (-1)**.*'^. 

eeC eeC eeC 

This proves the first point. To check the second one, observe that 
(_l)(«D-<dy)(a) = (-ifD(C)+l D ,{C)_ 

Clearly, £ D {C) + £ D '(C) = C ■ A{D, D') (mod 2), giving the proposition. □ 
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Corollary 4.3. Any dimer configuration D on a surface graph T C £ induces an 
isomorphism of affine -ff 1 (£;Z2)- spaces 

(3C(r)/~) -^Q(#i(£;Z 2 ),-), [if] g£ 
from the set of equivalence classes of Kasteleyn orientations onTcS onto the set 
of quadratic forms on (i?i(£; Z2), •)• Furthermore, ipn — ipD> if and only if D and 
D 1 are equivalent dimer configurations. 

Proof. The first part of Proposition 14.21 exactly states that -0£> is an isomorphism 
of affine Z2)-spaces. By the second part, iprj — ipD' if and only if the ho- 

momorphism -ffi(£; Z2) — ► 1i given by the intersection with A(D, D') is zero. By 
Poincare duality, this is the case if and only if A(D, D') — 0. □ 

4.2. Spin structures on surfaces. Recall that the fundamental group of SO(n) 
is infinite cyclic if n = 2 and cyclic of order 2 if n > 3. Hence, SO(n) admits a 
canonical 2-fold cover, denoted by Spin(n) — > SO(n). 

Let M be an oriented n-dimensional Riemannian manifold, and let Pso M 
be the principal SO(n)-bxmd\e associated to its tangent bundle. A spin structure 
on M is a principal Spin(n)-bwadle P — > M together with a 2-fold covering map 
P — > Pso which restricts to the covering map Spin(n) — > SO(n) on each fiber. 
Equivalently, a spin structure on M is a cohomology class £ E if 1 (-Pso! ^2) whose 
restriction to each fiber F gives the generator of the cyclic group H 1 (F; Z 2 ). 

It is well-known that such a spin structure exists if and only if the second Stiefel- 
Whitney class of M vanishes. In such a case, the set §(M) of spin structures on M 
is endowed with a natural structure of affine H 1 (M; Z 2 )-space. 

The 2-dimensional case is particularly easy to deal with for several reasons. First 
of all, any compact orientable surface £ admits a spin structure, as its second Stiefel- 
Whitney class is always zero. Furthermore, spin structures can be constructed using 
a certain class a vector fields, that we now describe. 

Let / be a non- vanishing vector field on £ \ a, where a is some finite subset of 
£. Recall that the index of the singularity x E a of / is defined as the degree of the 
circle map t 1— » /(7x(*))/|/(7x(*))|; where j x : S 1 — » £ \ a is a (counter-clockwise) 
parametrization of a simple closed curve separating x from the other singularities 
of /. Let us denote by V e „(£) the set of vector fields on £ with only even index 
singularities. We claim that any such vector field / defines a spin structure £/ on 
£. Indeed, consider a 1-cycle c in Pso (that is, a closed framed curve in £) and 
let us assume that c avoids the singularities of /. Then, let £/(c) E Z2 be the 
winding number modulo 2 of / along c with respect to the framing of c. Since all 
the singularities of / have even index, £/(c) = if c is a 1-boundary, and £/(c) = 1 
if c is a small simple closed curve with tangential framing. Therefore, it induces 
a well-defined cohomology class £/ E Hom(iJi (Pso ', Z2); Z2) = H 1 (Pso','^2) which 
restricts to the generator of the cohomology of the fibers. So £/ is a spin structure 
on £, and we have a map 

V e „(E) — »S(E), 

We shall need one last result about spin structures on surfaces, due to D. Johnson 
0. Given a spin structure £ E §(£), let q^ : -ffi(£; Z 2 ) — > Z 2 be the function defined 
as follows. Represent a E -ffi(£; Z 2 ) by a collection of disjoint regular simple closed 
curves 71,... , 7 m : S 1 c — > £. For all i and all i E S , complete the unit tangent 
vector ii{t)/\ii(t)\ to a positive orthonormal basis of T 7 .( t )£. This gives disjoint 
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framed closed curves in E, that is, a 1-cycle c in Pso- Set q^(a) = £(c) + m. 
Johnson's theorem asserts that is a well-defined quadratic form on (_ffi(E; Z2), •), 
where • denotes the intersection form. Furthermore, the map 

8(E)— ►Q(ffi(£;Z a ),.) ) 

is an isomorphism of affine _ff 1 (E; Z2)-spaces. 

4.3. Proof of Theorem 14.11 Let T C E be a surface graph, with E right-hand 
oriented. Given a Kasteleyn orientation K and a dimer configuration D on T C E, 
Kuperberg |TT| constructs a vector field /(if, D) £ V e „(£) as follows. Around each 
vertex of T, make the vectors point to the vertex. At the middle of each edge, make 
the vector point 90 degrees clockwise relative to the orientation K of the edge. 
Extend this continuously to the whole edges, as described in Figure El 




Figure 6. Kupcrberg's construction. 



The Kasteleyn condition J2J ensures that the vector field extends to the faces 
with one singularity of even index in the interior of each face. However, this vector 
field f(K) has an odd index singularity at each vertex of T. This is where the dimer 
configuration D enters the game: contract the odd index singularities in pairs along 
the dimers of D. The resulting vector field f(K, D) has even index singularities: 
one in the interior of each face of E, and one in the middle of each dimer of D. 

Gathering the results of the previous section and Kuperberg's construction, we 
get the following composition of maps: 

3C(r) x D(r) — V e „(E) — > 8(E) — Q(fr 1 (E;Z 2 ),-) , (K,D) q if(KD) . 

We are left with the proof that, given any Kasteleyn orientation K and dimer con- 
figuration D, the resulting quadratic form q = q^ J(K D) coincides with the function 
q^ defined in the statement of Theorem 14.11 So, given a G fli(E;Z2), represent 
it by a collection of simple closed curves C\, . .., C' m in T. (This is always possi- 
ble as r induces a cellular decomposition of E.) Since q is a quadratic form and 

« = E£ 1 [Ci], 

rn 

= (_i)E«<jC«-af]J(_i)«([cy). 
1=1 

Therefore, we just need to check that if C is an oriented simple closed curve in T, 
then (-l)'?([ c '])+ 1 = e K (C)(-1) £d ( c ). Consider the oriented regular curve 7 in E 
which follows C slightly on its left, and goes around the middle of each dimer it 
meets, except if it meets the same dimer twice. In this case, 7 stays close to C, as 
illustrated in Figure Clearly, 7 is a regular oriented simple closed curve in E. 
Furthermore, it is homologous to C and it avoids all the singularities of f(K,D). 
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We now have to check that the winding number u of f{K, D) along 7 with respect 
to its tangential framing satisfies (—1)" = e K {C){— 1)^ D ( C ). One easily check that 
to is equal (mod 2) to the winding number too of f(K, D) along 70, where 70 is the 
regular curve which goes around the middle of each dimer it meets, including the 
ones it meets twice. By construction of f(K, D), loq is equal to the winding number 
Co of f{K) along 7, where 7 is the regular curve which avoids all the dimers of D 
and all the vertices of T, as described in Figured 




Figure 7. The oriented simple closed curve C C T and its associ- 
ated regular curves 7, 70 and 7. The curve C is in solid, the dimers 
are in traced lines. 



The latter winding number Co can be computed locally by cutting 7 into pieces: 
one piece 7 e for each edge e of C, and one piece 7„ for each vertex v of T adjacent 
to a dimer of D sticking out to the left of C. The theorem now follows from the 
following case study. 

1. Let us assume that e is an edge of C such that £p(e) = +1. In this case, 
the vector field f(K) along 7 e in the tangential framing of 7 e defines a curve which 
is homotopically trivial, as illustrated in Figure |S] Hence, its contribution to Co is 
null. 

2. Consider now the case of an edge e of C such that s^(e) = —1. This time, 
f(K) along 7 e defines a simple close curve around the origin (see Figure IHJ). Its 
contribution to Co is equal to 1 (mod 2). 

3. Let v be a vertex of C with a dimer of D sticking out of v to the left of 
C. Then, the vector field f(K) along 7^ induces a simple closed curve around the 
origin, so its contribution to Co is equal to 1 (mod 2). The case illutrated in Figure 
121 is when K orients the dimer from v to its other boundary vertex. The other case 
is similar. 

Gathering all the pieces, the winding number of f(K) along 7 is equal to 
ILec £ c (e)(-lY D{c) - This concludes the proof of Theorem IO □ 

Using Johnson's theorem, we have the following immediate consequence of Corol- 
lary [0| 

Corollary 4.4. Any dimer configuration D on a surface graph T C £ induces an 
isomorphism of affine H 1 ^;^)- spaces 



i> D : (3C(T)/ ~) — S(£) 
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Figure 8. Computation of the local winding numbers. 



from the set of equivalence classes of Kasteleyn orientations on T C S onto the 
set of spin structures on E. Furthermore, ipu — ipD' if and only if D and D' are 
equivalent dimer configurations. □ 

5. PFFAFIAN FORMULAE FOR THE PARTITION FUNCTION AND CORRELATION 

FUNCTIONS 

5.1. The Kasteleyn matrix and its Pfaffian. Let if be a Kasteleyn orientation 
on a surface graph r C S with an even number of vertices. Enumerate these vertices 
by 1, 2, . . . , N = 2n. The Kasteleyn matrix is the 2n x 2n skew-symmetric matrix 
A K (T; w) = A K whose entry is the total weight of all edges from i to j minus 
the total weight of all edges from j to i. More formally, 

a£ = 5>tf( e M e )> 

e 

where the sum is on all edges e in T between the vertices i and j, e^(e) = if i = j, 
and 




1 if e is oriented by K from i to j; 
— 1 otherwise, 



if i ^ j. 

Consider a dimer configuration D on T given by edges e\, . . . , e n matching ver- 
tices ii and je for t = 1, . . . , n. It determines an equivalence class of permutations 
a: (l,...,2n) i— > . . . , i n ,jn) with respect to permutations of pairs (it, ft) 

and transpositions (ii,jt) i— > (je,ie). We will write this as a 6 D. Given such a 
permutation a, define 

s K (D) = (-iyf[s« je (e e ), 

where (— l) a denotes the sign of the permutation a. Note that this expression does 
not depend on the choice of a G D, but only on the dimer configuration D. 

Theorem 5.1. Let T CS be a surface graph with an even number of vertices. For 
any dimer configuration Dq on T and any Kasteleyn orientation K onT c£, 

(5) e K (D )Pf(A K )= Yl (-l) 9 °° {a) Z a (D Q ), 
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where q^ o is the quadratic form on -Hi(S; Z2) associated to K and Do, and 

D 

the sum being on all dimer configurations D such that A(Do,D) = a. 

Proof. Recall that the Pfaffian of a skew-symmetric matrix A — (a^ ) of size In is 
given by 

Pi(A) = (-l) CT <Ml)cr(2) • ■ • a ff (2n-l)<7(2n)) 

Men 

where the sum is on the set IT of matchings of {1, ... , 2n}. Therefore, 

Pi(A K ) = ^2 ( _ 1) C ' a a(l)a(2) ' ' ' a a(2n-l)a(2n) 

Men 

= ( _1 ) <T £ ?(D<T(2) ( e lM e l) ■■■Yl £ f(2n-l)<r(2n) {^n)w{e n ) 

[cr]en ei e„ 

n 

= Y^ 1 ^ II £ *(2i-i)a(2i) (ei)w(ei) 

D t=l 

= J2 £K (D)w(D), 

D 

where the sum is on all dimer configurations D on T. Hence, 
e K (D )Pf(A K ) =Y,z K {Do)e K {D)w{D). 

D 

Let us denote by e\, . . . , e„ the edges of T occupied by dimers of D, and by e\, . . . , e° 
the edges occupied by dimers of Dq. Fix permutations a and r representing the 
dimer configurations D and Dq, respectively, and set v — t o a . By definition, 

n n 

s K (D )s K (D) = (-l) T n £ Vi)W e °) ' (-^U^-iM^e) 

1=1 t=i 

n 

= ( — 1)" Y\. £ v(<j(2l-l))v(a(2t)){ e °z) £ a(2l~l)<j(2l)( e i)- 
1=1 

Note that the permutation v depends on the choice of a 6 D and t € Do, but it 
always brings the perfect matching D to Dq. Moreover, one can choose represen- 
tatives a G D and t E Do such that i/ is the counter-clockwise rotation by one 
edge of every (Do, Decomposition cycle C%, . . . , C m . For this particular choice of 
representatives, we have 

m m 

e K (Do)e K {D) = (-1)1^(^(^+1)^^ = 

i=l i=l 

Here, we use the fact that the length of a permutation cycle is the length of the 
corresponding composition cycle, and that the length of each composition cycle is 
even. Recall the quadratic form of Theorem 14.11 Since the Cj's are disjoint 
(Do, -Decomposition cycles, Ci ■ Cj = and ^D (Ci) = for all i, j. Therefore, 

m 

H(-e K (C t )) = (-l)^ A ^ D ». 

i=l 



16 



DAVID CIMASONI AND NICOLAI RESHETIKHIN 



The theorem follows. □ 

5.2. The partition function. Given a dimer configuration Do, Theorem pro- 
vides 2 2g linear equations (one for each equivalence class of Kasteleyn orientation) 
with 2 2g unknowns (the functions Z a (Do)). We want to use these equations to 
express the dimer partition function Z = Z a (Do) in terms of Pfaffians. 

Let V x V — * Z2, {ol,0) <— > a ■ (3 be a non-degenerate bilinear form on a Z2-vector 
space V. Recall that the Arf invariant of a quadratic form q: V — > Z2 on (V, •) is 
given by 

Arf(<z) = i£(-irH 

1 1 aGV 

Lemma 5.2. Let g, q' be two quadratic forms on (V, ■). Then, 

Arf{q)Arf{q') = (-l)^ = (-1)*'< A >, 
where A € V satisfies (q + q')(ct) — A ■ a for all a £ V . 

Proof. First note that q + q' is a linear form on V. Since the bilinear form ■ is 
non-degenerate, there exists A £ V such that (q + q')(a) = a ■ A for all a £ V. 
Furthermore, (q+q')(A) = A- A = 0, so (-1) 9 ( A ) = (-1) 9 '( A ). Let us now compute 
the product of the Arf invariants: 

Ar%)Ar%') = * £ (_i)«(«)4VW = _L £ (.^M^C/J+a). 
' ' q,/3sv ' ' q,/3gv 

Using the equality 

q(a) + q>{(3 + A) = q(a) + q'{(3) + g'(A) + (3 ■ A 

= q(a) + q(/3) + q'(A) = q{a + /?) + a ■ (3 + q(A), 

we obtain 



C_i w(A) 

Ar%)Ar%') = £(_i)«C*+fl+"0 



Q,/3 

= L M^(E 1 + E(- 1 ) 9(tt+ " )+Q ' /5 

We are left with the proof that the latter sum is zero: 

= ]T (-l)^ = ]T(-l)^V°-<), 

7^0 a+/3=7 7^0 

where is the cardinality of the set 

N* = {{a, (3) £ V x V \a + p = 7 and a ■ (3 = i} 

for i = 0,1. Since 7 7^ and (V, •) is non-degenerate, there exists x £ V such that 
x ■ 7 = 1. Then, the map (a,/?) 1— > (a + + .t) induces a bijection iV? — > iV*. 
Hence re = nz for all 7 7^ 0, and the lemma is proved. □ 

Using this lemma, it is easy to check that if dim(V) = 2n, then there are exactly 
2 2 ™" 1 + 2™ -1 quadratic forms on (V, •) with Arf invariant 1 and — 2™ _1 forms 

with Arf invariant —1. 
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Theorem 5.3. The partition function of a dimer model on a surface graph V C £ 
for a closed surface of genus g is given by the formula 

(6) Z= ^J2 Ar ^ )s K (D )Pf(A K ), 

[K] 

where the sum is taken over all equivalence classes of Kasteleyn orientations. Each 
summand is defined for a Kasteleyn orientation but it depends only on its equiva- 
lence class. Furthermore, the sign Arf(q^ g )e K (Do) does not depend on Dq. 

Proof. Recall that there is a free, transitive action of Z2) on the set of equiv- 

alence classes of Kasteleyn orientations onTcS ( Theorem 13. 20 . Let us denote by 
Kef, the result of the action of ^ G ff 1 (£;Z2) on a Kasteleyn orientation K. Then, 
by the first part of Proposition ^. 21 

(_l)(?o +«flo)W = (_!)*(») = : Xa (4>). 

Here Xa's are characters of irreducible representations of H l (Y,\ Z2). 
By equation ©, 

s K *{D )¥i{A K *)= (-l) q *° )ia) Z a (D ) 

oefli(S;Z) 

K 4> fa\ 

for all 4> £ Multiplying these equations by (— l) 9 - D o {p> and taking the 

sum over all <f> £ ff x (S; Z), we get 

(f) <fi a. 

= E(- 1 ) <) (a)+q °° m E (*)Sa (A>) • 

a ^ 

Using the orthogonality formula Xa(^)X/3(0) — 2 25 (5 a/ 3, we obtain 

MA)) = ^(-l) 9 "o( a! ^(5o)Pf(^). 
Therefore, the partition function Z = ^ Q Z a (Do) is given by 

where 

«t** = e **(A>)^ £ = e**(A,)Arf(«£)- 

a£ffl(S;Z 2 ) 

This gives the formula 

Z = 1 £ Arf(q^)e^ (Uo)Pf(^). 

0GH1(E;Z 2 ) 

Since -ff *(£, Z2) acts transitively and freely on the equivalence classes of Kasteleyn 
orientations, equality © follows. 

If K and K' are equivalent Kasteleyn orientations, then q^ o = q^ o so Arf (gjj ) = 
Ar%|Q. On the other hand, e K {D a ) = (-1)»e k '(Dq) andPf(A A ') = (-iyPf(A K ,), 
where ji is the number of vertices of T around which the orientation was flipped. 
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Therefore, the summand Arf(q^ o )e K (Do)Pf(A K ) does not depend on the choice of 
the representative in the equivalence class [K] . 

Let us finally check that the sign Ari(q^ g )e K (Dq) does not depend on Dq. Let 
D be another dimer configuration on T. By Proposition 14.21 Lemma |5. 21 and the 
proof of Theorem 15. II 

A(q« o )e K (D )A(q«)e K (D) = (-1)^ A <^» Arf(g£jA(g£) = 1. 
This concludes the proof of the theorem. □ 

5.3. Local correlation functions. In order to express local correlation functions 
< a ei ■ ■ ■ cr ek > as combinations of Pfaffians, let us recall several facts of linear 
algebra. 

Let A = {flij) be a matrix of size 2n. Given an ordered subset I of the ordered 
set a — (1, . . . , 2n), let Ai denote the matrix obtained from A by removing the i th 
row and the i th column for all i 6 I. Also, let (— l)""' denote the signature of the 
permutation which sends a to the ordered set I(a\I). If A is skew-symmetric, then 
for all ordered set of indices / = ■ ■ ■ ,ik,jk), 

d k Ff(A) (I 



® a hjl ' ' ' d a ikjk 



i-iyWpi(Ai). 



Furthermore, if A is invertible, then Pi(A) ^ and 

(-l) CT «Pf(A 7 ) = (-l^Pf^Pf^" 1 )^). 

So, let ei, . . . , efc be edges of the graph T, and let ie,je be the two boundary 
vertices of eg for I = 1, . . . , k. For simplicity, we shall assume that P has no 
multiple edges. Finally, set I = (ii,ji, ■ ■ ■ , ik,jk)- Applying the identities above to 
the Kasteleyn matrices A^ = A K *, Theorem 15.31 gives 



hT,' Kt U.<j,(- i r' ,) nAt 



If the Kasteleyn matrix is invertible for any Kasteleyn orientation of P, this 
expression is equal to 

k ^ 

^-xy* n<^ pf (^) pf ((^);\i)- 

(j> i=\ 

Since 

w(ei) ■ ■ -w{e k ) d k Z 



Z dw(e±) ■ ■ ■ dw(ek) ' 

the correlation functions are given by 

. fc E^*pf(^)n^^pf((^)-y 

<a ei ..-a ek >=(-!) ^—^—^ . 

if the Kasteleyn matrix is invertible for all possible Kasteleyn orientations of P. 
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If r is a planar graph, then the Kasteleyn matrix is always invertible since its 
Pfaffian is equal to the partition function. Therefore, 

< a ei ■ ■ ■ a ek >= (-l) k af in ■ ■ ■ a* ,Pf((^)^), 

where K is any Kasteleyn orientation on T C S 2 . 

On the other hand, there are graphs where some Kasteleyn matrix is not in- 
vertible. For example, consider a square lattice on a torus. Then, the Kasteleyn 
matrix corresponding to the spin structure with Arf invariant —1 is not invertible 
(see [IS]). 

Appendix A. Dimers and Grassman integrals 

A.l. Grassman integrals and Pfaffians. Let V be an n-dimensional vector 
space. Its exterior algebra AV = ©5J =0 A k V is called the Grassman algebra of 
V. The choice of a linear basis in V induces an isomorphism between AV and the 
algebra generated by elements fa, . . . ,<f> n with defining relations fa fa; = —fa; fa. 
The isomorphism identifies the linear basis in V with the generators fa. 

Choose an orientation on V. Together with the basis in V, this defines a basis 
in the top exterior power of V. The integral over the Grassman algebra of V of an 
element a £ AV is the coordinate of a in the top exterior power of V with respect 
to this basis. It is denoted by 

J a d(f). 

In physics, elements <fi are called Fermionic fields. More precisely, they are called 
neutral fermionic fields (or neutral fermions). They are called charged fermions if 
there is an action of U(l) on V. 

Recall that the Pfaffian of a skew symmetric matrix A of even size n is given by 

Pf ( A ) = ?¥ )n\ } (-ir^UMa) ' ' ' Mn-l)ff(n). 

where the sum is over all permutations a 6 S n . (This formula is easily seen to be 
equivalent to the one stated in Section 0) Expanding the exponent into a power 
series, one gets the following identity in AV: 

1 n 

Tr^exp (- ^2 OiAi A $7')) = Pf (^) 01 A • • • A <p n , 

where it: A V — > A n V is the projection to the top exterior power. In terms of 
Grassmann integral, this can be expressed as 

(7) / exp (~ o,o,,o,)do = Pf(A). 

Let us now assume that the space V is polarized, i.e. that V = W © W* with W 
some vector space and W* its dual. Then, the Grassman algebra of V is isomorphic 
to the tensor product of Grassman algebras for W and for W* in the category of 
super-vector spaces. That is, if the dimension of W is k, the Grassman algebra of V 
is isomorphic to the algebra generated by tpi, ijj*, i = 1, • • • k with defining relations 
ipi?pj = —Tpjiii, ipt>Pj — an d i>X^*j = —ipjipi- Such an isomorphism is 

specified by the choice of a linear basis in W. Note that such a choice induces a 
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basis in V (take the dual basis in W* ) and an orientation on V given by the ordering 
ipi, . . . , ip n , ip*, . . . , ip*. The Grassmann integral of an element a G Ay with respect 
to this choice of basis in W is denoted by 

adipdip* . 

Expanding the exponent, we obtain 

(8) / ex P ( E = (-l) M ^ Li Det(yl). 

Comparing this formula with (JJJ) , we obtain the following well known identity 
Pf f _ A t o J = (-l)^Det(A). 

Now consider the space U = V © V* . Let fa,...,(p n be a basis in V and 
fa , . . . , fa\ be the dual basis in V*. Using the change of variables Xi = til^EIij.^ 

X* = — y^"*^ - , together with the equalities (JTJ and ijSJ, we obtain 



exp (- a H^i + 2 E o-tiVi^. 
= (-i) -1 ^ y ex P ( E XiatjX*j)dxdx 



= Det(A). 

One can easily derive other Pfaffian idcntitcs in a similar way . 

A. 2. Dimer models on graphs and Grassman integrals. Let r C E be a 

surface graph, and let A K = [a^ ) be the Kasteleyn matrix associated with a 
Kasteleyn orientation K on T C E. By Theorem 15.31 and identity (0), the partition 
function for dimers on T is given by 



Z = h E Arf«)e K (^o) | exp (1 £ ^ 



[K] '' '~t,jev(r) 



where V(r) denotes the set of vertices of T. 

Recall that e K (D ) = (— 1) CT n"=i £ f e j e > where the dimer configuration D matches 
vertices ig and ji for t = 1, . . . , n, and a is the permutation (1, . . . , In) i— > (ii, ji, . . . , i n ,jn)- 
Taking into account the identity 

n n 

Mr n £ %jt #i • • ■ = n dfa e dfa e , 

f=l £=1 

the formula for the partition function can be written as 

z = ^E[^(l E 

[K] J i,jev(r) 

where 

n 

D K (f> = Atf(q« Q )e K (D ) dcP = Ar%£ ) J] e*. dfa e dfa e . 
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Similarly, local correlation functions can be written as 
1 1 ^ 

< <y ei •••<r efc >= — / s &°«&) n^^^^ 

[X] J 1=1 

Here and j; are the boundary vertices of the edge e;, and 7 = . . . , ik,jk)- 
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